Local phonon mode in a fermionic bath, and its relation to the Kondo effect 
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We have studied the interplay of a local phonon mode embedded in a metallic host (Holstein 
impurity model) using Abelian bosonization. The phonon frequency softens, which takes place 
in two steps: first, their frequency starts softening, and acquires finite lifetime. Then oscillations 
disappear from the response, and two distinct, finite dampings characterize them. Similar behaviour 
shows up in the spin-boson model. Due to phonons, the electrons experience an attractive, dynamic 
interaction. As a result, the electronic charge response enhances similarly to the spin response in 
the Kondo model. The local electronic density of states develops a dip around zero frequency. Thus 
the chance of charge-Kondo effect emerges. 



PACS numbers: 73.23-b,73.63-b,72.10.Fk 

I. INTRODUCTION 

Single impurity models have a long historjsi, and our 
understanding has profited a lot from different reliable 
techniques. Numerical renormalization group, conformal 
field theory, Nozieres' Fermi liquid description, Bethe 
ansatz solution and bosonization completed each other, 
and highlighted different aspects of the same problem. In 
many cases, the latter, being nonperturbative and ana- 
lytic, was able to provide us with a transparent picture of 
the underlying physics^. Its appeal is due the elegance of 
analytic solution, being able to treat dynamic quantities 
as well (correlation functions), and the simplicity of the 
picture emerging from it. 

Beyond electron-electron interaction, which is a basic 
ingredient in Kondo type models, another source of cor- 
relation is represented by phonons^iii^i^iii^. They play 
a prominent role in explaining the conventional s-wave 
BCS superconductivity and the charge density wave for- 
mation in low dimensional systems. Their general feature 
is the generated dynamic interaction between the elec- 
trons, which, in essence, is attractive for small energy 
transfers. 

In many strongly correlated systems, such as heavy 
fermions or valence fluctuation systems, lattice vibrations 
are known to couple strongly to electrons. The recent dis- 
covery of magnetically robust heavy fermion behaviour 
in filled skutterudite compound^ SmOs4Sbi2 renewed in- 
terest in Kondo phenomena with phononic originiSili. 
Both the specific heat coefficient and coefficient of the 
quadratic temperature dependence of the electrical resis- 
tivity were found to be almost independent of an applied 
magnetic field. Quantum dots, especially the ones based 
on single molecules (Cgo for examplei^), also possess vi- 
brational degrees of freedom, which will react to the elec- 
tron transfer through them^ Therefore the properties of 
impurity models dominated by phonons are challenging. 
The applied methods should be able to treat dynami- 
cal properties, which rules out techniques as the Bethe 
ansatz. 

Several theoretical investigation focused on electron 
transport in the presence of electron-phonon interaction. 



mainly in the context of the single impurity Anderson- 
Holstein model. This model consists of an Anderson 
impurity with linear coupling to a local phonon mode^. 
Hewson and Newns considered its spinless version with 
a few electron in the system without explicit electron 
correlatio n^^' Later further studies based on the NRG 
(numerical renormalization group)^'^'^^ and the non- 
equilibrium Keldysh formalisni^i^ have been performed. 
These works reported about the softening of the local 
phonon mode, and the enhancement of the charge suscep- 
tibility. These phenomena point toward the realization 
of the charge-Kondo effect, caused by the degeneracy of 
zero and doubly occupied electron states. 

In the present work, using Abelian bosonization, for 
the first time to our knowledge to attack the local 
electron-phonon problem, we can not only confirm the 
prediction of previous works, but also study analyti- 
cally the specific heat, the phonon Green's function, the 
charge susceptibility and especially the local density of 
states. We follow the softening of phonons from weak to 
strong coupling, and connect the present problem with 
the underscreened Kondo model in magnetic field. As 
the electron-phonon coupling increases, phonon states 
with large occupation number start to play an impor- 
tant role even at very low temperatures, leading to a 
transition to a phonon distorted state only at infinitely 
strong coupling. Similar phenomenon reveals itself in the 
spin-boson model as well^^'^^. Then, in contrast to NRG, 
where only a finite number of phonons can be taken into 
account, we can consider states with arbitrary phonons, 
and the crossover and the softening at arbitrary electron- 
phonon couplings. 



II. THE MODEL 

In real materials, electrons interact with each other 
explicitly through the Coulomb force and implicitly 
through phonons. In several cases, the resulting 
phase is well described by Landau's Fermi liquid the- 
ory with renormalized quasiparticle parameters due to 
interactionsiS. However, certain inhomogeneities like im- 
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purities, lattice imperfections are usually present, and ex- 
hibit a source of scattering. When these are of electronic 
origin, spin-Kondo physics is evoked^. The local per- 
turbations can also show up in the form of local lattice 
vibrations. This is why we have undertaken the study 
of spinless fcrmions interacting with a single, dispersion- 
less Einstein phonon mode, mainly at zero temperature. 
This can be called the single impurity Holstein model, 
which is related to the Holstein model as the single im- 
purity Anderson model is related to its periodic version. 
It also parallels to the one studied in Ref. fl^lTsI in the 
same way as the single impurity Wolff model^S is con- 
nected with the Anderson model^. A similar model has 
been studied^^ describing the coupling of a local Einstein 
phonon to conduction electrons by the numerical renor- 
malization group, with emphasis on the nature of the 
fixed points, and relations with the two-level Kondo ef- 
fect have been revealed. The present model can also be 
relevant for the dynamical mean-field theory (DMFT) 
of the Holstein lattice model^i, which maps it onto the 
Holstein impurity model. Coupling to acoustic phonons 
would also be a fruitful proposition, but the essence of 
physics is readily captured by the simplest model. When 
the coupling of the electron density to the phonon dis- 
placement field is isotropic, the model can be mapped 
onto a single branch of chiral fermions interacting with a 
single Einstein phonon only at the origin, and is suitable 
for Abelian bosonization. This mapping follows closely 
the one proposed in Ref. for the Kondo model-. 
In the Hamiltonian language, the model is given by: 
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and only the radial motion of the particles is accounted 
for by chiral (right moving) fermion field'^^, p{x) —: 
"i'^{x)'^(x) : is the normal ordered electron density, L 
is the length of the system, v is the Fermi velocity, g de- 
scribes the local electron-phonon coupling, m and luq are 
the phononic mass and frequency, respectively, Q and P 
are the phonon displacement field and momentum con- 
jugate to it. 

In spirit, this model is similar to the underscreened 
Kondo model in a magnetic field (which quenches the re- 
maining degrees of freedom) with phonons replacing the 
impurity spin, since the phonon displacement field (Q) 
can take any real values, while the electron density try- 
ing to compensate it, is bounded, although states with a 
large number of bosons in the harmonic oscillator hardly 
contribute to the physics at low temperatures due to their 
high energy. The underscreened Kondo model in a mag- 
netic field is governed by a Fermi liquid fixed point. We 
speculate that the model under study produces similar 
behaviour. When the magnetic field is switched off, the 
underscreened Kondo model shows singular Fermi liq- 
uid behaviour—, whose analogue in the present case oc- 
curs at g — > cxD. Since the charge degrees of freedom 



are coupled to the local bosons, the possibility of ob- 
serving the " underscreened" charge-Kondo effect opens^. 
The realization of the charge-Kondo effect has also been 
found in similar model o^^'^^ . Our model is similar to the 
one describing a tunneling particle coupled to a fermionic 
The model can be bosonized via32i2^ 
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If we integrated out the phonon degrees of freedom in Eq. 
(H]) or ([3]) , we would arrive to a local interaction between 
electrons, given by 2g'^ /m{uj'^ ~ ojq), uj is the energy 
transfer in the interaction. At low energies (|a;| < luq), 
this would lead to an effective attractive interaction, and 
the resulting Hamiltonian would coincide with the spin 
sector of the repulsive Wolff impurity mode P°'^^i'^° , re- 
sponsible for the spin-Kondo phenomenon. Hence our 
model is capable to show Kondo physics in the charge 
sector. This further strengthens our proposal for the 
charge-Kondo effect. Attractive local interactions has al- 
ready provided us with similar phenomenon^'^ . 

This resulting Hamiltonian (Eq. ([3])) is similar to the 
Caldeira-Leggett Hamiltonian (CL), describing the dy- 
namics of a particle coupled to dissipative environment'^"'^ . 
However, the roles are reversed between oscillators and 
tunneling particle in Eq. ([3]) and in Ref. [3l|: our local- 
ized phonon represents the tunneling particle in CL lan- 
guage, while the bosonized (<&) fermionic field stands for 
the dissipative environment of CL, represented by har- 
monic oscillators. An important difference between the 
two models is the lack of renormalization of potential in 
CL, ensuring, that the coupling to reservoir solely intro- 
duces dissipation^^. In our case, however, we expect the 
softening of the phonon mode due the interaction with 
fermions on physical ground, hence the renormalization 
of the phonon frequency plays an essential role as is tes- 
tified later in Eq. (|14p . in addition to dissipation. More- 
over, the properties of our original fermionic field must 
be obtained through Eq. a difficulty which is never 
faced with in the Caldeira-Leggett model. 

Our bosonized Hamiltonian (Eq. ([3])) also shares sim- 
ilar properties with the spin-boson model, where a two- 
level system is coupled to a bosonic environmen t ^ "^i . 
There, assuming an Ohmic bath, at a finite value of 
the coupling of the boson-two-level system, a quantum 
phase transition has been identified. In our case, the 
bosonic bath is represented by the bosonized fermions, 
but the two-level system is replaced by a multi-level sys- 
tem, i.e. a harmonic oscillator. As a result, the ob- 
tained behaviour of our model differs from that found in 
the spin-boson model, although some dynamical features, 
connected with the oscillation frequency, are similar. 
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III. LIMITING CASES 

The adiabatic limit (luq ^ 0) is meaningful, if it is 
taken together with m — > cx), keeping mujQ = k finite. 
The phonons kinetic energy vanishes, and only a static 
displacement is felt by the electrons. Hence, the model 
reduces to a simple scattering problem, although one has 
to minimize the total free energy of the system (elec- 
tronic contribution + elastic energy with respect 
to the frozen phonon displacement field to obtain the cor- 
rect solution, which transforms our problem into a self- 
consistent mean-field theory, which is the exact solution 
in this limit. The free energy, subject to minimization 
with respect to Q, reads at T = as 

1 

2 f K 

AF = -nogQ— / dX a.Tcta.n{\TT gQ p) + -Q'^ , (4) 
IT J 2 


which leads to the self-consistency equation for Q: 

Q — arctaniTTgQp), (5) 
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where uq — Wp is the average electron density per 
site. In the weak coupling limit, only the trivial so- 
lution exists {Q — 0), as opposed to standard mean- 
field treatments of the Holstein lattice model^. As g 
exceeds gc = yKj2nQp, a non-trivial solution shows 
up, which minimizes the free energy. This behaves as 

, 3/2 

Q — \J'o{g~ gcijirpgc close to gc-, and Q = gno/n as 
g grows to infinity. For g > gc, the electrons feel a lo- 
cal scattering center with strength V = gQ, and the local 
electronic density of states is expected to be suppressed22, 
close to zero frequency as 1/(1 + {nVp)'^). Complete sup- 
pression occurs only for g — > oo. Phonons lose their dy- 
namics completely. However, when quitting the strict 
adiabatic limit, the displacement field is not frozen any 
more, quantum corrections are expected to destroy or- 
dering. In essence, this is similar to the large N limit of 
impurity models^ (N denotes the spin degeneracy), where 
a true phase transition occurs only in the 1/N=0 limit. 
Note, that similar calculation applies to the adiabatic 
limit of the Anderson-Holstein impurity model as well. 

In the anti- adiabatic limit, luq —^ oo while keeping 
the ratio g/ujo finite, the phonons react instantaneously 
to the electrons, leading to non-retarded interaction be- 
tween the densities as 
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where n is the charge density of electrons at the ori- 
gin. Had we chosen spinful fermions, the resulting model 
would coincide with the negative U Wolff modcP"-'^'^ 
{U = —g'^/mcuQ). In the spinless case, this again re- 
duces to potential scattering with strength V = g'^{2no — 
l)/2mLLjQ. Again, the local density of states is expected 
to be suppressed close to zero frequency, similarly to the 
adiabatic case. 



Finally by taking the atomic or zero bandwidth limit, 
which also describes the infinite coupling case {g — > oo), 
the phonons couple only to an isolated electron: 



Hat = c+c{E - 
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where E is the c-level energy. This model is known as the 
independent boson modelS. Since the electron number 
operator (c^c) is conserved, most of the physical quan- 
tities can trivially be evaluated. However, the electron 
Green's function does not belong to this class. After a 
unitary transformation {U ~ e'x:p{~igPc^c/mujQ)), this 
can be evaluated exactly as 
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Here the first term is the renormalized c-level Green's 
function, (Q^) = l/2mujQ is the mean square value of the 
displacement field, being independent of the coupling in 
this case. Finally, DQ{t) is the Green's function of the 
displacement field Q given by 



DQ{t) = (g2) cxp(-zc^oi). 



(9) 



From this, one can evaluate the density of states, which 
reads as 



Pat{ 
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with 7 = g"^ /2muj^. Compared to the original density of 
states with g = 0, where a single Dirac-delta peak con- 
tains all the spectral weight, it is now distributed non- 
uniformly among infinite number of peaks. In general, 
most of the spectral weight is concentrated around uj k E 
for a finite g as well. In the weak coupling limit (7 <C 1), 
the weight decreases monotonically with n. For strong 
couphng (7 » 1), the largest weights are still distributed 
around lo E, which means that the n « 7 terms are 
the most dominant in the sum. These weights are pro- 
portional to \l ^pn. From this, one can conclude, that by 
increasing g to move from weak to strong coupling, terms 
with n cx 7 contain the largest spectral weights, which, 
however decreases as 1/ ^^7. This means that increasing 
amount of spectral weight is transfered away from the 
central region, hence the local density of states is sup- 
pressed due to phonons, as borns out from the above 
three limiting cases. It is only suppressed to zero for 
infinitely strong coupling g. 



IV. SOFTENING OF THE PHONON 

To determine the dynamics of our system, we start 
with the evaluation of the Green's function of the phonon 
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displacement field, defined by 

Dgir) = (r,Q(T)Q(0)), (11) 

whose Matsubara form can easily be calculated using the 
standard diagrammatic technique from Eq. ([3]), which 
yields to 

1 1 

Dqiiujm) = — 2 — o 2 t (12) 

TO ^4 2pg^x{^m)/m 

where p = 1 /27rw, ujm is the bosonic Matsubara frequency 
and 

xM = E 2 ^rr ^2 (13) 

is the local susceptibility of the fermions without the 
phonons. After analytic continuation to real frequencies, 
it is given by x(ijJ ^ W) « p{W + iTruj/2) for small fre- 
quencies, W is the cutoff or bandwidth, and the obtained 
formula holds regardless to the chosen cutoff procedure. 
Our Eq. (|12p is exact for the Hamiltonian given by Eq. 
([3]), no other diagrams are left to consider. However, 
the bosonization cutoff scheme renormalizes the coupling 
constant^, and the impurity potential g needs to be re- 
placed by the appropriate scattering phase shift. Since 
this is not a universal quantity (apart from the linear 
term), it depends on the regularization scheme, and we 
consider it to be replaced by the renormalized one on- 
wards, as is standard in impurity models^. The dynamics 
of the phonons can be inferred by investigating the pole 
structure of Eq. The excitations energies become 

complex (i.e. finite lifetime or damping of phonon excita- 
tions) in the presence of finite coupling to the electrons, 
and are given by 

ujp± = -ir ± y'w^ - r2 - ATW/n (14) 

with r = 7r(5p)V2m. Below E < Ei -2T4^/7r + 
^/4W^J^T^ + uJ^, the square root in Eq. (fT4|) is real, hence 
the displacements reach their equilibrium in an oscilla- 
tory fashion within a time characterized by — 1/ImWp. 
For higher frequencies, the phonon mode is completely 
softened, ReiVp = (the square root becomes imagi- 
nary), the excitations have two different finite lifetimes 
or dampings determined from Eq. (|14p between Fi < 
r < r2 = ttujq/AW. Interestingly, the very same criti- 
cal value r2 would be found, if we inserted the value of 
= clIo\/ m/2nop found in the extreme adiabatic limit 
(or self-consistent mean field theory) to F = 7r(gp)2/2m, 
which is, in turn, zero because loq — > 0. In this range, 
phonons can be excited with zero energy, and all the 
displacements are relaxed to equilibrium without oscilla- 
tions. This relaxation slows down close to F2. This re- 
gion is very narrow, because for realistic values, loq <^ W, 
it shrinks as ~ uiq/W^. Similar relaxation characterizes 
two-level systems coupled to a bosonic batbiiii^. For 
higher values of F, our approach breaks down, it signals 



lattice distortion with (Q) ^ 0, as is indicated by the 
complete softening of the phonon mode. However, such 
a phenomenon is impossible in zero dimensional systems, 
and we ascribe it to bosonization. The explicit functional 
form of physical quantities on model parameters obtained 
by bosonization can deviate from the exact one. The 
phase shifts in impurity problems'^^, the correlation ex- 
ponents in Luttinger liquids determined via bosonization 
are only correct in the weak coupling limit^. This and 
the previous argumentation suggests, that in reality, F2 
should only be reached for g —^ 00, but the complete soft- 
ening predicted at Fi would take place at a finite value 
of g, and F ~ only at small g. In the adiabatic limit, 
(to — > 00, luq ~f 0) T also vanishes, and lu± — 0. 

We mention that such softening remains absent in the 
Caldeira-Leggett model due to the compensation of the 
potential renormalization thcrc'^^. Similar phenomenon 
occurs in a two-level system, coupled to a bath of three 
dimensional phonons"^ The renormalized tunneling rate 
drops to zero only at infinitely strong coupling. 

The softening of the phonon frequency was also found 
in the NRG treatment of the Anderson-Holstein impurity 
model^i^ii^, but this has only been followed in a moder- 
ate range of parameters. DMFT studies of the Holstein 
lattice model also revealed similar features^i. 

The general behaviour of the eigenfrequencies is 
sketched in Fig. [TJ The phonon density of states contains 
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FIG. 1: (Color online) The real (blue solid line) and imaginary 
(red dashed line) part of the phonon excitation energies, and 
the inverse of the mean square value of the displacement field 
(black dashed-dotted line) are visualized as a function of the 
electron-phonon coupling for W/uq = 2. At Pi, the phonon 
modes is softened, but the excitations have finite lifetime. At 
the damping vanishes, the lifetime (~ — l/Imojp) diverges, 
and so does the mean square value of the displacement field. 

two Lorentzians with resonance widths -Imwp-t , centered 
around ReWp^. With increasing F, they become centered 
around the origin, and close to F2 , the one determined by 
LUp-f- tends to 6{lu). This is reminiscent to what happens 
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in the channel anisotropic two channel Kondo model^^. 
In comparison with the underscreened Kondo model in 
a magnetic field, a similar divergence shows up in the 
localized electron density of states^'^ with the offset of 
magnetic field. The very same role is played in our case 
by r — > which causes the divergence as ujp+ 0. 
Softening of the phonon mode has been found by the nu- 
merical renormalization group as weli^ in a similar model. 



V. THERMODYNAMICS 

The mean value of the displacement field is zero. The 
mean square value of the displacement field is calculated 
from Eq. ^ at T = as 



2iujQ 



■In"^- 



(15) 

where Qq — l/2mujQ is the value for g = 0. It is plotted 
in Fig. [H As r approaches it diverges as 



Ql 



4ri 



(16) 



indicating the change in the (Q) — relation. This is 
traced back to the transition to the phonon distorted 
state at (7 = cx), as discussed below Eq. pij) . Similar phe- 
nomenon has been observed in the spin-boson modeU^ii^ 
at a finite value of the spin-boson coupling, assuming 
an Ohmic bath. However, for a super-Ohmic bath, the 
transition occurs only at infinitely strong coupling^^. In 
contrast, in Caldeira-Leggett-type Hamiltonians, dissipa- 
tion induced squeezing occurs'^^, which originates from 
the effect of damping in reducing the rate of escape from 
a metastable well. 

In a linear harmonic oscillator, the ground state wave 
function is a gaussian. It is natural to ask, to what extent 
this picture holds in the case of finite electron-phonon 
coupling. The probability distribution of the oscillator 
coordinator or displacement reads as 



\^oscix)\^ - {S{Q-X)) 



lim 

a—fO 



dk 
2^ 



(17) 

In the second step, a useful representation of the Dirac- 
delta function was inserted. Since the Hamiltonian is 
quadratic, or the corresponding action is gaussian for the 
phonons after integrating out the fermions, the expecta- 
tion value of the exponent can be calculated following 
Ref. m, which leads to 



|*osc(2:)|' 



1 



: exp 
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The ground state wave function remains gaussian, but 
the variance ((Q^)) increases monotonically with g. The 
region where the phonons are mainly restricted to, is 



wider than without electron-phonon coupling. The cou- 
pling to electron increases the average oscillator displace- 
ment around the equilibrium. This further corroborates 
the picture emerging from the previous studies. 

The phonon contribution to the free energy can be cal- 
culated by using Pauli's trick of integrating over the cou- 
pling constant. After some algebra, one arrives at 



dx 



2rx 



2n 



b(x) tan ""^ , „ „ , 

^' ' x2-c.2(i-(r/r2)) 
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where fig is the free energy of the fermions in the ab- 
sence of phonons, b{x) is the Bose distribution function, 
and the integral should be limited to the bandwidth, but 
this can be sent to infinity when obtaining quantities by 
differentiating fi"^^. One has to carefully chose the ap- 
propriate phase angle of the tan~^ function. This ex- 
pression is similar to the free energy of the two chan- 
nel Kondo model in a magnetic field (~ lvq \/i~-(r/T2) 
here) along the Emery-Kivelson line^^, after replacing the 
Bose distribution function with the Fermi one. How- 
ever, our "magnetic field" cannot be switched off, hence 
the entropy is zero at T = 0, similarly to the under- 
screened Kondo modeP^ . The specific heat can be calcu- 
lated by C(T) = -l/T{d'^n/dT^). Its low temperature 
(T ^ ujQy^l — (r/r2)) behaviour due to phonons reads 
as 



Cp{T) = 



ttT r 
6^r2-r' 



(20) 



which sharpens when r2 is approached. The exponential 
freezing-out of the Einstein phonon changes to a linear 
T dependence, and adds to the Sommerfeld coefficient of 
the conduction electrons, and enhances it. Had we cho- 
sen acoustic phonons, their specific heat would also 
be overwhelmed at low temperatures. At higher temper- 
atures, a broad bump shows up in the specific heat as a 
function of temperature due to the presence of g. This 
can serve as an identifier of the local electron-phonon in- 
teraction, and is shown in Fig. (5] Similar structures have 
been observed in the Andcrson-Holstein impurity model 
by NRGii. 



VI. ELECTRONIC PROPERTIES: 
SUSCEPTIBILITY AND LOCAL DENSITY OF 
STATES 

The phonons are expected to have a profound impact 
on the electronic properties. Had we integrated them 
out, we would have arrived to a local dynamic electron- 
electron interaction, attractive at low energies and repul- 
sive for higher ones. Static interactions already modify 
significantly the electronic response, as was demonstrated 
in the Wolff model2£i2£ii2,. In the foUowings we are going 
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FIG. 2: (Color online) The total (blue solid line) phonon spe- 
cific heat contribution for different values of the coupling con- 
stant r/r2 = 0.1, 0.3, 0.5, 0.7 and 0.9, VF/wo = 10 from 
bottom to top. The red dashed dotted line shows the pure 
phononic F = specific heat, while the black dashed curves 
stands for their difference due to the finite coupling g. Note 
the sharpening of the broad bump at low T as F increases. 



to see the effect of dynamic one explicitly on the electron 
system. 

The electric charge is defined through p{x) = 
dx^{x) I "v/tt. The local charge correlation function is eval- 
uated from Eq. ([3]) using the diagrammatic technique, 
and along the real frequency axis it reads as 



Xchargeiy') 



1 



(21) 



x(w) 



The local charge excitation energies are determined from 
the poles Eq. (HH), which possess finite lifetime, as given 
by Eq. (fT4)l . Eq. ([2T|) is the standard result in the 
random-phase-approximation, which turns out to be the 
exact one in the present case. Most of the many body 
contributions are canceled by the Ward identity, which 
relates the vertex function and the electron propagator. 
This could be expected from the fact that the model is 
solvable by bosonization. The very same phenomenon 
revealed itself during the study of the Wolff impurity 
model22ii^. The dynamic nature of the phonons is ob- 
servable in the denominator of Eq. (|2ip . In the static 
limit, the effective interaction between the electrons is 
attractive, and changes to repulsive for frequencies ex- 
ceeding the phonon energy ujq. The static limit of the 
charge susceptibility simplifies to 



Xc/iar(/e(0) — 



2p^W 



(i-(r/r2))- 



(22) 



In accordance with Eq. (IMI) , this also predicts the transi- 
tion to the distorted phase, which would be accompanied 



by the rearrangement of the charges as well. This signals 
that our model is on the brink of charge Kondo effect. 
Similar enhancement of the charge response has been re- 
ported in Ref. [HI 

Since we work with spinless electrons, the single 
particle Green's function can be evaluated, unlike in 
the SU(N) Wolff model^Si^S or the two channel Kondo 
model^S. There, such a calculation would involve for- 
mally yJW(x), which is difficult to work with. The local 
retarded Green's function is defined as 



G«(^) = -^e(^)({M'(^),M'+(0)}) 



exp(-47r($2)) 
27ra 



[exp(47rC(t)) 



-exp(47rC(-i))], 
(23) 



where in the second step we made use of the gaussian na- 
ture of the action and used the usual tricks of operator 
manipulation^^'—. The task is to evaluate the correla- 
tor C{t) — ('i>(t)$(0)), whose second appearance in Eq. 
([23]) follows from time reversal symmetry. Given the fact 
that the bosonized Hamiltonian (Eq. ([3])) is quadratic, 
we can evaluate this expectation value at bosonic Mat- 
subara frequencies using the path integral representation 
of the problem-S. Then by transforming it to imaginary 
times, finally we can read off C{t) after careful analytic 
continuation to real times. First we find for the Matsub- 
ara form that 



1 



-TDQ[Uyn). 



(24) 



Here the first term is responsible for the 1/r decay of 
fermionic correlations. Fortunately, the second expres- 
sion, accounting for the phonon contribution, is separated 
from the first one, and does not require ultraviolet reg- 
ularization unlike the first one, where the large lo„i part 
need to be cut off to avoid divergences. 

Compared to the exact expression derived in the 
atomic hmit {Gat{t) in Eq. ([SD), Eqs. d^Sll-lEll) have 
identical structures. It contains the free fermionic con- 
tribution, and in addition, the exponentiated phonon 
Green's function. Since the coupling constant appears 
in the exponent, together with Dq, similarly to the ex- 
act solution in the atomic limit, this suggests the non- 
perturbative nature of our bosonization approach in the 
presence of a finite electron band. Hence we believe, that 
our study captures correctly the physics in the whole pa- 
rameter range. 

Its Fourier transform with respect to a;,„ yields to 
the imaginary time ordered expression for the correlator, 
from which the desired relation can be obtained as 



C(i)-C(0) = -^lnf^- 
47r \ a + ivt 



-T{DQ{t)-{Q^)l 
(25) 

where the equal time correlator need to be subtracted to 
regularize the first term on the right hand side and 

2i 



Dqit) 



mniLUp^ — ojp-) 



{f{uj.p+t) - fiiUp^t)) (26) 
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with 

7r 

f{x) = — sin(x)Si(x) — — sin(a;) — cos(x)Ci(— x), (27) 

where Si(a;) and Ci(a;) are the sine and cosine integrals. 
At large times, its decay is characterized by the damp- 
ing, lmLUp± , and the frequency of oscillations by ReWp± . 
By plugging this formula to GR{t), we evaluate the re- 
tarded single particle Green's function. The local density 
of states follows as 

oo 

p{ui) — Im / dtex.p(iu!t)Gii{t),— 

TT J 

(oo 
l + zexp(-2.™r(Q^))|^H^P|^x 


X [exp(27rmri:)Q(-t)) - exp{2TrniTDQ{t))]) . (28) 

It is shown in Fig. [3] for various values of F. As F 
increases from weak (F/r2 <C 1) to strong couphng 
(r/r2 ^ l)j the residual density of states decreases 
monotonically. In addition, increasing amount of spec- 
tral weight is transfered away from the central region 
around w = 0, in perfect agreement with what we found 
in limiting cases of the model in Sec. IIIII Complete sup- 
pression at = occurs in all cases at infinitely strong 
coupling, similarly to the bosonized case at F — > F2 . This 
is why we believe that our approach captures correctly 
the strong coupling limit as well. Close to the critical cou- 
pling F2 , the density of states exhibits a V-shapes form 
around the Fermi energy. Similar redistribution of spec- 
tral weight was identified in the adiabatic, anti-adiabatic 
and atomic limit. These findings are also in accordance 
with those found by NRG in the Anderson- Holstein im- 
purity model^i^i^. There, in the non-interacting version 
of the model, a resonance peak shows up at the f-level 
energy, which narrows significantly in the presence of 
phonons, and increasing amount of spectral weight is 
transfered away from this region. A more direct com- 
parison is difficult, since the Anderson- Holstein impu- 
rity model distinguishes between conduction and local- 
ized electrons, which are weakly hybridized with the con- 
duction band, as opposed to our one-band model. Sim- 
ilar reasoning separates the purely electronic Anderson 
impurity model from the Wolff impurity model'^*^. 

The significant decrease (dip hke structure) in p(0) 
should be detected in point contact spectroscopy or by 
scanning tunneling microscopy experiments in bulk ma- 
terials. For small values of F, a steplike drop occurs 
with decreasing frequency close to uj — tjQ-^/l — F/F2, 
which smoothens as the electron-phonon coupling in- 
creases, through the appearance of other steplike fea- 
tures, as is seen in Fig. [3]for F — O.7F2. As ujq increases, 
the drop in the density of states occurs at higher fre- 
quencies (oc uJQ^y^^^T/T2)^ and the suppression of the 
residual density of states at w = is more efficient for 
higher wq, as can be seen in the inset of Fig. [3) 



These features are different from that induced by a 
non-magnetic impurity or by a magnetic one, and the 
reason is the effective, dynamically attractive interaction 
generated by the phonons. 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 



uj/oJo 

FIG. 3; The change of the local electron density of states at 
the impurity site is shown for W/loq = 10 from weak to strong 
coupling for r/r2=0.1 (black), 0.3 (red), 0.5 (blue dashed- 
dotted), 0.7 (blue dashed) and 0.9 (blue solid) from top to 
bottom. Increasing amount of spectral weight is transfered 
away from the zero frequency region with P. The inset visual- 
izes the change in the residual density of states for W/luq = 10 
(solid) and 20 (dashed), which vanishes at F = 



VII. CONCLUSION 

In conclusion, we have studied a model of spinless 
fermions interacting with an Einstein phonon (single im- 
purity Holstein model) using Abelian bosonization. With 
increasing coupling, the phonon mode softens, and at 
g — > 00, distortion occurs. The softening of the phonon 
resonance takes place in two steps: first, the frequency 
starts softening, and acquires finite lifetime. Then os- 
cillations disappear from the response, and two distinct, 
finite dampings characterize the phonons. Finally, at a 
critical coupling, which is conjectured to be at g = 00, the 
phonons are distorted. These manifest themselves in the 
specific heat through the enhancement of the Sommerfeld 
coefficient. The phonon's entropy is similar to that of the 
Kondo model in finite magnetic field'^^, hence vanishes at 
T = 0. The present model resembles closely to the un- 
derscreened Kondo model in magnetic field. The local 
charge susceptibility is strongly enhanced, similarly to 
the spin susceptibility of the Kondo model. This points 
toward the realization of charge-Kondo effect. Closed 
expression is derived for the local electronic density of 
states. With increasing coupling, significant amount of 
spectral weight is transfered away from the uj = region. 
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in accordance with results in the atomic hmit and in the 
Anderson-Holstein impurity model. This should be ob- 
servable in point contact spectroscopy, or on a molecule 
trapped near a tunnel junction, and influence the current- 
voltage characteristic as well. 
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